Strategies in crowd and crowd structure 
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In an emergency situation, imitation of strategies of neighbours can lead to an order-disorder phase 
transition, where spatial clusters of pedestrians adopt the same strategy. We assume that there are 
two strategies, cooperating and competitive, which correspond to a smaller or larger desired velocity. 
The results of our simulations within the Social Force Model indicate that the ordered phase can 
be detected as an increase of spatial order of positions of the pedestrians in the crowd. 
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I. INTRODUCTION 

Simulations of crowd dynamics is among the most ac- 
tively developed applications of physical methods in so- 
cial sciences. Standard difficulties met in such applica- 
tions are in this case either absent, or seriously weak- 
ened. First, crowd dynamics is relevant for control of 
large gatherings and prevention of stampede catastro- 
phes [l], 0. Second, experimental data are accessible 
[1,13. Third, a model (Social Force Model, SFM) has 
been developed [BJ which can be used by physicists with- 
out any input of social sciences. Fourth, realistic values 
of the model parameters have been obtained and verified 
Q. These causes made SFM a popular tool, used and 
developed by various scientific groups of physicists and 
computer scientists. To develop its connections with col- 
lective sociological effects is one of our motivations here. 

In the original version of SFM , most of the param- 
eters are related to physical aspects of the interactions be- 
tween pedestrians and of the interactions of pedestrians 
with obstacles. There are two exceptions: the strength of 
interpersonal, psychologically motivated repulsion A and 
the desired velocity Vd- A variation of A has been used 
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in 0, Q to include individual decisions on calling for help 
and coming to aid. Here we are interested in a variation 
of the desired velocity Vd- In SFM, equations of motion 
for each pedestrian are solved numerically, and one of the 
forces is proportional to the difference Vd — v(t), where 
v(t) is the actual velocity of the given pedestrian. In the 
conditions of evacuation, i>d=1.8 m/s seems to be appro- 
priate as an accepted velocity of brisk walking [{J. On 
the contrary, Wd=3.0 m/s was found to produce jamming 
at narrow exits [HJ. Here we use these two values to 
differentiate between two strategies of evacuation: coop- 
erating (1.8 m/s) and competitive (3.0 m/s). 

The collective character of the process is due to the 
imitation effect. Namely, we use probabilities w(n) that 
the strategy (cooperative or competitive, X = C or S) 
of a given pedestrian will be changed if n out of 6 di- 
rect neighbours of a given pedestrian adopt this strategy. 
The number six is the number of nearest neighbours in 
a triangular two-dimensional lattice, what is a model of 
people in crowd; suggested by our former observations 
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[lfj. In our former work, we used a square lattice 
yet we do not expect qualitative differences between the 
results for these two structures. 

Our aim here is as follows. We expect that the imita- 
tion leads to spatial correlations of strategies of pedestri- 
ans, and that these correlations can lead to a respective 
phase transition [ll|. We ask if this phase transition is 
visible in the crowd structure. This interest is motivated 
by the idea of detecting collective effects online by anal- 
ysis of crowd motion, registered by local cameras. 

The strategy of our modelling is as follows. First we set 
the probabilities w(n) as dependent on a single parameter 
x and we solve numerically the Ising-likc problem of se- 
lecting strategies C or S by pedestrians. As the outcome 
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we get the set w c (n) where there is a transition between 
the ordered and the disordered phase. Next we verify 
the obtained transition point by means of the fluctua- 
tion function [12]. These steps are described in Section 
HTl Next we perform the simulations of crowd dynam- 
ics within SFM, where pedestrians select Vd= 1.8 m/s 
(strategy C) or 3.0 m/s (strategy S). During the simu- 
lation, we evaluate how a local crowd structure depends 
on the parameter x and, in particular, we check if the 
above mentioned phase transition is visible. This step is 
described in Section IIIII These considerations are illus- 
trated by the results of our field observations of crowd, 
described in Section IIV1 Last section is devoted to dis- 
cussion. 



II. THE PHASE TRANSITION 

At this preliminary stage of modelling pedestrians are 
placed at nodes of a two-dimensional, triangular lattice 
and their motion is not taken into account. The lat- 
tice contains N = 10 6 sites and helical boundary con- 
ditions are assumed. In the initial state, strategies C 
and S are assigned to pedestrians with given probabili- 
ties r, 1 — r; spatial correlations are not taken into ac- 
count. During the simulation, pedestrians change their 
strategies with some probabilities w(n), that depend on 
the number n of neighbouringsubjects currently using 
the same strategy (see also Q). Here the functions 
w(n) depend on a single control parameter x as follows: 
w(0) = 1, w(l) = 3x, w(2) = 2x, w(3) = x, w(4) = x/2, 
w(5) = x/4, w(6) = x/6. If w{n) > 1 for any n, we set 
w(n) = 1. For simplicity, these probabilities are the same 
for both strategies C and S. 

The numerical results indicate that indeed an order- 
disorder phase transition appears near x c = 0.429, as 
shown in Fig. [T] For x < x c , the probabilities of state 
changes are small and clusters of pedestrians with the 
same strategy persist. This result is confirmed by the 
analysis of temporal dependence of the order parame- 
ter, defined as m = (N(C) - N(S))/N; there, N(X) 
is the number of agents who adopt strategy X, and 
N = N(C) + N(S) = 10 6 . The temporal dependence 
of m for several control parameters x = 0.45, 0.5 and 0.6 
(x > x c ) is presented on Fig. [2] The simulation takes 
M = 10 8 steps and single step is completed when all 
N = 10 6 nodes arc investigated in typewriter order. The 
last half of these steps (M/2 = 5 ■ 10 7 ) is used for average 
order parameter (m) evaluation (Fig. [!}. 

The continuous character of the phase transition is con- 
firmed by studying the persistence of order parameter 
time series m(t) using multifractal detrended fluctation 
analysis with linear detrending (MF-DFA1) 12j. In this 
approach, the q-ih order fluctuation function is defined 
as 
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where F 2 (is,s) is the detrended variance calculated for 
each segment v of the length s, which is obtained by the 
division of the analysed time series into non-overlapping 
segments of equal length; q can be any real value, here 
we use q e {-2, 0.5, 2}. 
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FIG. 1: Average order parameter (m) dependence on control 
parameter x. The last M/2 = 5 • 10 7 steps of m(t) temporal 
evolution is used for average (m) evaluation. 
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FIG. 2: (Colour online) The temporal dependence of order 
parameter m(t) for several values of the control parameter 
x — 0.45, 0.5 and 0.6 (x > x c ). The range of fluctuations 
increase when we approach the transition point x c from above. 



Figure [3] shows typical MF-DFA fluctuation functions 
for x w 0.42. One can see a crossover from random walk 
type scaling behaviour (F q (s) ~ s hl - q } with h(q) « 1.5) 
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on small time scales s to random white noise type scaling 
behaviour (h(q) « 0.5) on large time scales. The tempo- 
ral range T of the correlations is given by the position of 
the crossover, being close to 10 4 in this case. Figure 0] 
shows the dependence of T on x, which follows a power 
law, T(x) ~ (x - Xc) 13 with [3 « -2. We see that T con- 
tinously diverges in the vicinity of the transition point. 
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FIG. 3: MF-DFA fluctuation functions F q (s) versus the scale 
s for q £ {—2, 0.5, 2}; as it is seen, the curves for different q's 
almost overlap. 




To inspect the crowd structure, for each pedestrian i we 
identify his six nearest neighbours and we calculate the 
difference Di between the distances from i to his nearest 
neighbour and from i to his sixth nearest neighbour. In 
the triangular lattice, the difference is zero. The smaller 
Di, the crowd is more ordered. 

In Fig. [51 we show the dependence of the mean differ- 
ence of distances (D) on the parameter x, measured at 
the tentatively stationary state. The mean is calculated 
over all pedestrians i = 1, • • ■ ,N and over a time period 
of about ten seconds, after a transient time of the same 
length. As we see, (D) is larger but more stable above 
the transition point, i.e. for x > x c . Below this value 
of £, a clear reduction of (D) can be observed, although 
large fluctuations hamper the accuracy. On the other 
hand, the information on the initial state, encoded by 
the probability r, is preserved below the transition; this 
is shown in Fig. [5] 
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FIG. 5: The mean difference of distances {D) as dependent 
on x, taken from the SFM simulations. 



FIG. 4: The temporal range T of correlations as dependent 
on x. 



III. THE CROWD SIMULATIONS 

The simulation is performed in accordance to the orig- 
inal formulation of the Social Force Model and the model 
parameters Pedestrians are moving, all in the same 

direction, along a corridor of five meters width, with pe- 
riodic boundary conditions. After each one second, they 
check the strategies adapted at the moment by their six 
nearest neighbours, and they modify their strategies ac- 
cording to the probabilities w(n), the same as given in 
Section HU 



IV. FIELD OBSERVATION 

The Corpus Christi procession in Krakow, Grodzka 
Street, was registered by a hand camera on June 7, 2012 
from a window at 1-st floor, about five meters above the 
street. The movie lasts ten minutes, but after about first 
five minutes and 30 seconds the crowd density started to 
decrease. Also, during two first minutes priests and mil- 
itary orchestra arc moving in a prescribed order. Then 
the material is limited to about 3 '30", when the crowd 
density is about one person per m 2 . 

The aim was to check if some structure of the crowd 
can be found. For this purpose, a lane along the street 
is selected of about three meters width. For pedestri- 
ans who appear at this lane, the times when their heads 
disappear from the screen are registered manually. The 
differences between these times are used to make a his- 
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FIG. 6: The mean value (r) of the percentage of the pedestri- 
ans who adopt the competitive strategy, taken from the SFM 
simulations. Below x c , the results depend on the initial value 
of r, which was taken as 0.7 (upper curve) and 0.3 (lower 
curve) . 



togram of time gaps r, shown in Fig. [7] We note that 
the velocity of the crowd is more or less constant. 
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FIG. 7: The experimental histogram of time gaps r between 
pedestrians, as taken from the movie. 



We expected and hoped to see two peaks in the his- 
togram of the time gaps r, and indeed the obtained plot 
does not shatter these expectations. First peak, at about 
0.2 s, is seen in the plots for all movies we did, and it can 
be safely interpreted as related to pedestrians who walk 
together. The second peak about 0.5 s is specific for this 
crowd density. Once the density decreases, the times be- 
tween subsequent groups vary much more and the second 
peak is not visible. 

V. DISCUSSION 

The results of numerical simulations within SFM indi- 
cate, that if the imitation process leads to a phase transi- 
tion, and if the two states are related with different values 
of the desired velocity, then the ordered phase manifests 
itself by an increase of short-range order. On the con- 
trary, in the disordered phase the spatial correlations de- 
crease, what means that clusters of people moving with 
the same velocity disappear. 

It would be desirable to confirm this conclusion by our 
field observation. However, the data recorded are re- 
lated only to an ordered phase, where pedestrians move 
with approximately the same velocity, and the statistics 
is rather poor. Then, the obtained spectrum of times be- 
tween neighbouring pedestrians can serve here as merely 
an illustration that some amount of spatial short range 
order does exist. 

Our results suggest, that it is possible to infer about 
the state of crowd from its spatial structure registered 
in one instant of time. More local ordering indicates the 
presence of collective effects, which can be an effect of 
imitation. This tool of crowd control can be complemen- 
tary to correlations of velocity of pedestrians, detected 
from movies 
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